A host-microbiota interactome reveals extensive transkingdom connectivity.

The myriad microorganisms that live in close association with humans have diverse effects on physiology, yet the molecular bases for these impacts remain mostly unknown1-3. Classical pathogens often invade host tissues and modulate immune responses through interactions with human extracellular and secreted proteins (the 'exoproteome'). Commensal microorganisms may also facilitate niche colonization and shape host biology by engaging host exoproteins; however, direct exoproteome-microbiota interactions remain largely unexplored. Here we developed and validated a novel technology, BASEHIT, that enables proteome-scale assessment of human exoproteome-microbiome interactions. Using BASEHIT, we interrogated more than 1.7 million potential interactions between 519 human-associated bacterial strains from diverse phylogenies and tissues of origin and 3,324 human exoproteins. The resulting interactome revealed an extensive network of transkingdom connectivity consisting of thousands of previously undescribed host-microorganism interactions involving 383 strains and 651 host proteins. Specific binding patterns within this network implied underlying biological logic; for example, conspecific strains exhibited shared exoprotein-binding patterns, and individual tissue isolates uniquely bound tissue-specific exoproteins. Furthermore, we observed dozens of unique and often strain-specific interactions with potential roles in niche colonization, tissue remodelling and immunomodulation, and found that strains with differing host interaction profiles had divergent interactions with host cells in vitro and effects on the host immune system in vivo. Overall, these studies expose a previously unexplored landscape of molecular-level host-microbiota interactions that may underlie causal effects of indigenous microorganisms on human health and disease.

The gut microbiome modifies the associations of short- and long-term physical activity with body weight changes.

BACKGROUND: The gut microbiome regulates host energy balance and adiposity-related metabolic consequences, but it remains unknown how the gut microbiome modulates body weight response to physical activity (PA). METHODS: Nested in the Health Professionals Follow-up Study, a subcohort of 307 healthy men (mean[SD] age, 70[4] years) provided stool and blood samples in 2012-2013. Data from cohort long-term follow-ups and from the accelerometer, doubly labeled water, and plasma biomarker measurements during the time of stool collection were used to assess long-term and short-term associations of PA with adiposity. The gut microbiome was profiled by shotgun metagenomics and metatranscriptomics. A subcohort of 209 healthy women from the Nurses' Health Study II was used for validation. RESULTS: The microbial species Alistipes putredinis was found to modify the association between PA and body weight. Specifically, in individuals with higher abundance of A. putredinis, each 15-MET-hour/week increment in long-term PA was associated with 2.26 kg (95% CI, 1.53-2.98 kg) less weight gain from age 21 to the time of stool collection, whereas those with lower abundance of A. putredinis only had 1.01 kg (95% CI, 0.41-1.61 kg) less weight gain (pinteraction = 0.019). Consistent modification associated with A. putredinis was observed for short-term PA in relation to BMI, fat mass%, plasma HbA1c, and 6-month weight change. This modification effect might be partly attributable to four metabolic pathways encoded by A. putredinis, including folate transformation, fatty acid ß-oxidation, gluconeogenesis, and stearate biosynthesis. CONCLUSIONS: A greater abundance of A. putredinis may strengthen the beneficial association of PA with body weight change, suggesting the potential of gut microbial intervention to improve the efficacy of PA in body weight management. Video Abstract.

High-sensitivity pattern discovery in large, paired multiomic datasets.

MOTIVATION: Modern biological screens yield enormous numbers of measurements, and identifying and interpreting statistically significant associations among features are essential. In experiments featuring multiple high-dimensional datasets collected from the same set of samples, it is useful to identify groups of associated features between the datasets in a way that provides high statistical power and false discovery rate (FDR) control. RESULTS: Here, we present a novel hierarchical framework, HAllA (Hierarchical All-against-All association testing), for structured association discovery between paired high-dimensional datasets. HAllA efficiently integrates hierarchical hypothesis testing with FDR correction to reveal significant linear and non-linear block-wise relationships among continuous and/or categorical data. We optimized and evaluated HAllA using heterogeneous synthetic datasets of known association structure, where HAllA outperformed all-against-all and other block-testing approaches across a range of common similarity measures. We then applied HAllA to a series of real-world multiomics datasets, revealing new associations between gene expression and host immune activity, the microbiome and host transcriptome, metabolomic profiling and human health phenotypes. AVAILABILITY AND IMPLEMENTATION: An open-source implementation of HAllA is freely available at http://huttenhower.sph.harvard.edu/halla along with documentation, demo datasets and a user group. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Strain Identification and Quantitative Analysis in Microbial Communities.

Microbiology has long studied the ways in which subtle genetic differences between closely related microbial strains can have profound impacts on their phenotypes and those of their surrounding environments and communities. Despite the growth in high-throughput microbial community profiling, however, such strain-level differences remain challenging to detect. Once detected, few quantitative approaches have been well-validated for associating strain variants from microbial communities with phenotypes of interest, such as medication usage, treatment efficacy, host environment, or health. First, the term "strain" itself is not used consistently when defining a highly-resolved taxonomic or genomic unit from within a microbial community. Second, computational methods for identifying such strains directly from shotgun metagenomics are difficult, with several possible reference- and assembly-based approaches available, each with different sensitivity/specificity tradeoffs. Finally, statistical challenges exist in using any of the resulting strain profiles for downstream analyses, which can include strain tracking, phylogenetic analysis, or genetic association studies. We provide an in depth discussion of recently available computational tools that can be applied for this task, as well as statistical models and gaps in performing and interpreting any of these three main types of studies using strain-resolved shotgun metagenomic profiling of microbial communities.

Bayesian modelling of high-throughput sequencing assays with malacoda.

NGS studies have uncovered an ever-growing catalog of human variation while leaving an enormous gap between observed variation and experimental characterization of variant function. High-throughput screens powered by NGS have greatly increased the rate of variant functionalization, but the development of comprehensive statistical methods to analyze screen data has lagged. In the massively parallel reporter assay (MPRA), short barcodes are counted by sequencing DNA libraries transfected into cells and the cell's output RNA in order to simultaneously measure the shifts in transcription induced by thousands of genetic variants. These counts present many statistical challenges, including overdispersion, depth dependence, and uncertain DNA concentrations. So far, the statistical methods used have been rudimentary, employing transformations on count level data and disregarding experimental and technical structure while failing to quantify uncertainty in the statistical model. We have developed an extensive framework for the analysis of NGS functionalization screens available as an R package called malacoda (available from github.com/andrewGhazi/malacoda). Our software implements a probabilistic, fully Bayesian model of screen data. The model uses the negative binomial distribution with gamma priors to model sequencing counts while accounting for effects from input library preparation and sequencing depth. The method leverages the high-throughput nature of the assay to estimate the priors empirically. External annotations such as ENCODE data or DeepSea predictions can also be incorporated to obtain more informative priors-a transformative capability for data integration. The package also includes quality control and utility functions, including automated barcode counting and visualization methods. To validate our method, we analyzed several datasets using malacoda and alternative MPRA analysis methods. These data include experiments from the literature, simulated assays, and primary MPRA data. We also used luciferase assays to experimentally validate several hits from our primary data, as well as variants for which the various methods disagree and variants detectable only with the aid of external annotations.

Functionalization of CD36 cardiovascular disease and expression associated variants by interdisciplinary high throughput analysis.

CD36 is a platelet membrane glycoprotein whose engagement with oxidized low-density lipoprotein (oxLDL) results in platelet activation. The CD36 gene has been associated with platelet count, platelet volume, as well as lipid levels and CVD risk by genome-wide association studies. Platelet CD36 expression levels have been shown to be associated with both the platelet oxLDL response and an elevated risk of thrombo-embolism. Several genomic variants have been identified as associated with platelet CD36 levels, however none have been conclusively demonstrated to be causative. We screened 81 expression quantitative trait loci (eQTL) single nucleotide polymorphisms (SNPs) associated with platelet CD36 expression by a Massively Parallel Reporter Assay (MPRA) and analyzed the results with a novel Bayesian statistical method. Ten eQTLs located 13kb to 55kb upstream of the CD36 transcriptional start site of transcript ENST00000309881 and 49kb to 92kb upstream of transcript ENST00000447544, demonstrated significant transcription shifts between their minor and major allele in the MPRA assay. Of these, rs2366739 and rs1194196, separated by only 20bp, were confirmed by luciferase assay to alter transcriptional regulation. In addition, electromobility shift assays demonstrated differential DNA:protein complex formation between the two alleles of this locus. Furthermore, deletion of the genomic locus by CRISPR/Cas9 in K562 and Meg-01 cells results in upregulation of CD36 transcription. These data indicate that we have identified a variant that regulates expression of CD36, which in turn affects platelet function. To assess the clinical relevance of our findings we used the PhenoScanner tool, which aggregates large scale GWAS findings; the results reinforce the clinical relevance of our variants and the utility of the MPRA assay. The study demonstrates a generalizable paradigm for functional testing of genetic variants to inform mechanistic studies, support patient management and develop precision therapies.

Design tools for MPRA experiments.

Motivation: Genetic reporter assays are a convenient, relatively inexpensive method for studying the regulation of gene expression. Massively Parallel Reporter Assays (MPRA) are high-throughput functionalization assays that interrogate the transcriptional activity of many genetic variants at once using a library of synthetic barcoded constructs. Despite growing interest in this area, there are few computational tools to design and execute MPRA studies. Results: We designed an online web-tool and R package that allows for interactive MPRA experimental design encompassing both power analysis and design of constructs. Our tool is tuned using data from real MPRA studies. Users can adjust experimental parameters to examine the predicted effect on assay power as well as upload VCFs for automated construct sequence generation. Availability and implementation: The MPRA Design Tools web application is available here: https://andrewghazi.shinyapps.io/designmpra/, https://github.com/andrewGhazi/designMPRA and https://github.com/andrewGhazi/mpradesigntools. Supplementary information: Supplementary data are available at Bioinformatics online.